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Abstract 

We review various B meson decays that require knowledge of the transverse 
decay constant of the fei(1235) meson. We report on an exploratory lat- 
tice QCD calculation of the transverse decay constant of the b\ meson. The 
lattice QCD calculations used unquenched gauge configurations, at two lat- 
tice spacings, generated with two flavours of sea quarks. The twisted mass 
formalism is used. 
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1. Introduction and motivation 

The transverse decay constant of the &i(1235) meson (/^) is theoretical 
input to a number of decays of the B meson. For example is an important 

QCD input to the following decays: B° — > bip + , B — > bjK* + , using the 
light cone formalism [Tj [2]. The constant is also input to the decay 
B — > bij P H] that use light cone sum rules. Diehl and Hiller [5] discuss 
studying decays of the B meson with final states that include the b\ meson. 

There are alternative theoretical formalisms [HI U\ El E] to the light cone 
sum rules such as factorization, that describe the non-leptonic decays of the 
B meson to final states that include the bi meson. So it is important to have 
a cross-check of the input parameters used in the light cone formalism. 
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BaBar has experimentally measured the B decays: b\R and b\K [TUl [TTj . 
The charmless decays of the B meson, that include those with a b\ meson in 
the final state, have been reviewed by Cheng and Smith [T2] . 

The transverse decay constant of the bi meson is not accessible to exper- 
iment, but can be calculated in models [13] and sum rules. Calculations of 
the bi meson have also been used to tune sum rules [TJ1 [T3J [TB] . In particular, 
the same sum rules are used to simultaneously extract the transverse decay 
constants of the bi and p mesons [HI [15] . 

In principle lattice QCD should be able to produce an accurate result for 
particularly as modern lattice QCD calculations usually have multiple 
lattice spacings and volumes, with pion masses below 300 MeV [17]. To the 
best of our knowledge, there has never been a lattice QCD calculation of /? 
before this one. 

The b\ meson is a good test case for lattice techniques that deal with 
resonances, because it is thought to be a basic quark-antiquark meson that 
decays via S-wave. The experimental width of the &i(1235) is 142(9) MeV, 
and the bulk of the decays are to un. Hence further motivation for this study 
is to compute as much information about the b\ meson from our lattice QCD 
calculations as possible. 

Light cone sum rules and factorization methods, also use the decay con- 
stants of the do, 7r(1300), a\ mesons to study the decays of the B me- 
son, but there have been previous lattice QCD calculations of those quanti- 
ties psiCSlED!. 

2. The lattice QCD calculation 

The transverse decay constant (/^(//)) of the b\ meson is defined [21] by 

(0 | I h(P, A)) = i/J(/i)e^ef A) P^ (1) 

where a^ u = i/2[j fl , j v ], and e?U is the polarisation vector of the meson. It is 
convenient to introduce the tensor current T Ufl = ipa^ip. We do not include 
any momentum in this lattice calculation. For completeness we note that the 
&i meson has J PC = l + ~. 

In the isospin limit the leptonic decay constant of the b\ meson is zero, 
because the operator is orthogonal to the vector and axial currents. Al- 
though for some quantities, such as the matrix element for p — u mixing [22] 
or the decay constant of the flavour non-singlet ++ meson [19] an estimate of 
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Tabic 1: Summary of results for fjf used in this calculation. afi q is the bare mass of the 
light quark in lattice units. The ensemble names are from |28j . 

isospin violating quantities can be made from a lattice calculation with two 
degenerate flavours of sea quarks, we don't see how to estimate the leptonic 
decay constant of the b\ meson without using non-degenerate light quarks. A 
formalism to do this has recently been developed for twisted mass QCD [23J . 

Our lattice calculation uses the twisted mass QCD formalism [23]. Once 
a single parameter has been tuned, twisted mass QCD has non-perturbative 
0(a) improvement [25, 26J. The twisted mass formalism has been reviewed 
by Shindler [27]." 

We have recently reported on the some basic measurements of the p, b± 
and ao mesons [29J. In this paper we extend that study to the transverse 
decay constant of the b\ meson. All the necessary details are in the previ- 
ous paper [2S] and here we provide a brief summary. There is a additional 
information about the lattice techniques, such as the smearing and varia- 
tional analysis in the "methods paper for the ETM collaboration [30]." We 
used the twisted mass Wilson action and the tree level improved Symanzik 
action. The ensembles used in this analysis are in table [2] Details of the 
analysis of light pseudo-scalar mesons are in [201 EI] • 

The correlators used to extract the decay constant are in equation [2 
We use a smearing matrix of order 2, using basis functions of local and fuzzec 
interpolating operators, that includes the correlators in equation [2j We fit the 
smearing matrix to a factorizing fit form [30J with two states. The charged 
interpolating operator for the b\ meson in the twisted basis was used [50] • 
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Figure 1: Effective mass plot for the b\ channel at f3 = 3.9, fi q =0.004, L=24. The labels 
LL, FL, FF refer to local-local, fuzzed-local, and fuzzed- fuzzed. 
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E E e^TufoOfyMt)^— *^ (2) 

x k=l,i<j 

We used ensembles at two different (3 values. At f3 = 3.9 we included two 
volumes. We used the pion decay constant of the 7r meson to determine the 
lattice spacing. At f3 — 4.05 (/3 = 3.9) the lattice spacing is: a = 0.0667(5) 
fm (a = 0.0855(5) fm). The lattice spacing from f n was consistent with the 
value from the nucleon mass [32] ■ In figures [l] and [2] we show effective mass 
plots for the b\ correlators. 

The decay constant fjf depends on the value of the renormalisation scale. 
We used a renormalisation factor obtained from the Rome- Southampton non- 
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Figure 3: The decay constant of the b\ meson as a function of the square of the pion 
mass. 

perturbative method [33j EH EH] . As traditional in lattice QCD calculations 
we quote the result at the scale of 2 GeV. The renormalisation group equa- 
tions can be used to evolve the decay constant to another scale. 

The results for the decay constant from this calculation are in table [2] In 
figure [3] the transverse decay constant of the b\ meson is plotted in physical 
units as a function of the square of the pion mass. Although the error bars 
are large, the results for are consistent between the two lattice spacings 
and volumes. 

A common technique to check whether a state is a scattering state or a 
resonance is look at the volume dependence of the amplitude [35] ■ The use 
of the volume dependence of the amplitude and the connection to Luscher's 
method has recently been discussed by Meng et al. [37] . The "rule of thumb" 
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is that if our b\ correlator couples to a scattering state of urn, the volume 
dependence of the decay constant extracted from equation [2] is 

fl ~ ^ (3) 

where V is the spatial volume. For a resonance fjf should be independent 
of the volume (apart from small corrections if the box size is too small to fit 
the resonance state). 

The numerical results at \i q = 0.004 at /3 — 3.9, show that f^(L = 
32)//^ (L = 24) = 0.94(24), compared to the prediction for scattering states 
in equation [3] of 0.65. 

In our previous paper |29j we showed that the decay of the b\ meson to 
un was open in our calculation. However the mass of the lightest state in 
the 1 h channel didn't track the sum of the masses of the u and ir mesons 
particularly well. In Liischer's formalism for the study of resonances on 
the lattice the opening of strong decays is described by an avoided level 
crossing [38J . The detailed calculations of Bernard et al. [31] for the A baryon 
suggested that the avoided level crossing is "washed out" by the dynamics, 
so comparing the mass from the resonant interpolating operator to the sum 
of the masses of the decay products is probably too simplistic. A similar 
situation happened with string breaking, where the linearly rising potential 
was seen to increase beyond the energy that allowed the string to break to 
two static B mesons [4"0] . 

The above considerations suggest that although we don't have full control 
over the resonant nature of the 61 meson, our interpolating operators are 
coupling to the b\ meson in the range of quark masses in our calculations. 
We extrapolate fjf linearly in the square of the pion mass to get /^(2 GeV) 
= 236(23) MeV at the physical pion mass at j3 — 3.9. 

In table [2] we collect together other estimates for We have evolved 
our lattice result to the scale of 1 GeV to compare with the results from sum 
rules. The formalism to evolve the decay constant with scale is described [29J , 
that uses input from perturbative calculations by Gracey and others HH1 
ETJ EH]- The perturbative factor is 1.095 to evolve fT from 2 GeV to 1 GeV. 

In the non-relativistic quark model the b\ meson is a P-wave meson with 
a node in the wave-function at the origin. This would suggest the transverse 
decay constant is small. The partial inclusion of relativistic effects in the 
quark model [H] increases the decay constant. The results in table [2] show 
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Bakulev and Mikhailov 021 g3] 


sum rule 


181(5) 


Yang gl EQ 


sum rule 


180(8) 


This calculation 


lattice QCD 


258(25) 



Table 2: Summary of calculations of at 1 GeV. 



that f? is of the same order of magnitude as the pion decay constant (132 
MeV), so this is evidence that local interpolating operators will couple to 
the bi meson. Pragmatically using derivative sources [501 151~] with smearing 
techniques, such as Jacobi, may be useful to get a good signal. 



3. Conclusion 

We have presented the first calculation of the transverse decay constant 
of the 6i meson from lattice QCD. We obtain /£(2 GeV) = 236(23) MeV at 
the physical pion mass. The result is higher than the sum rule results by 3 
a. Future lattice QCD calculations need to reduce the statistical errors on 
the correlators, and to directly take into account the resonant nature of the 
h\ meson. 

4. Acknowledgments 

We thank all members of ETMC for a very fruitful collaboration. We 
thank Prof. Braun for providing the error for the result in [UJ, and we 
thank Mihail Chizhov for useful comments on the paper. 

References 

[1] K.-C. Yang, Phys. Rev. D72, 034009 (2005), [hep^ph/0506040| 



[2] H.-Y. Cheng and K.-C. Yang, Phys. Rev. D76, 114020 (2007), 
0709.0137, 



[3] M. Jamil As lam and Riazuddin, Phys. Rev. D75, 034004 (2007), |hep"- 
|ph/0607114| 



8 



[4] W. Wang, R.-H. Li, and C.-D. Lu, (2007), 0711.0432, 



[5] M. Diehl and G. Hiller, JHEP 06, 067 (2001), hep-ph/0105194 



[6] V. Lap orta, G. N ardulli, and T. N. Pham, Phys. Rev. D74, 054035 

(2006) , |hep^ph/0602243[ 

[7] G. Calderon, J. H. Munoz, and C. E. Vera, Phys. Rev. D76, 094019 

(2007) , 0705.1181, 

W. Wang, R.-H. Li, and C.-D. Lu, Phys. Rev. D78, 074009 (2008), 
0806.2510, 

[9] J. H. Munoz and N. Quintero, (2009), 0905.1526, 

[10] The BABAR, B. Aubert et al, Phys. Rev. Lett. 99, 241803 (2007), 
0707.4561, 

[11] BABAR, B. Aubert et al, Phys. Rev. D78, 011104 (2008), 0805.1217, 

[12] H.-Y. Cheng and J. G. Smith, (2009), 0901.4396, 

[13] M. V. Chizhov, JETP Lett. 80, 73 (2004), [hep^ph/0307100 , 



[14] A. P. B akulev and S. V. Mikhailov, Eur. Phys. J. C17, 129 (2000), 
|hep-ph/9908287| 

[15] W. Broniowski, M. V. Polyakov, H.-C. Kim, and K. Goeke, Phys. Lett. 



B438, 242 (1998), hep-ph/9805351 



[16] J. Govaerts, L. J. Reinders, F. de Viron, and J. Weyers, Nucl. Phys. 
B283, 706 (1987), 

[17] K. Jansen, (2008), 0810.5634, 

[18] M. Wingate, T. A. De Grand, S. Collin s, and U. M. Heller, Phys. Rev. 
Lett. 74, 4596 (1995), |hep^ph/9502274| 

[19] UKQCD, C. McN eile and C. Michael, Phys. Rev. D74, 014508 (2006), 
|hep-lat/0604009[ 

[20] UKQCD, C. Mc Neile and C. Michael, Phys. Lett. B642, 244 (2006), 
|hep-lat/0607032| 



[21] K.-C. Yang, Nucl. Phys. B776, 187 (2007), 0705.0692, 

[22] European Twisted Mass, C. McNeile, C. Michael, and C. Urbach, Phys. 
Lett. B674, 286 (2009), 0902.3897, 

[23] A. Walker-Loud, (2009), 0904.2404, 

[24] Alpha, R. Frezzotti , P. A. Grassi, S. Sint, and P. Weisz, JHEP 08, 058 
(2001), |hepTat/010100H 

[25] R. Frezzotti and G. C. Rossi, JHEP 08, 007 (2004), hep-lat/03060l"4] 

[26] P. Dimopoulos, R. Frezzotti, C. Michael, G. C. Rossi, and C. Urbach, 
(2009), 0908.0451, 

[27] A. Shindler, Phys. Rept. 461, 37 (2008), 0707.4093, 

[28] European Twisted Mass, C. Urbach, PoS LAT2007, 022 (2007), 
0710.1517, 

[29] K. Jansen, C. McNeile, C. Michael, C. Urbach, and f. t. E. Collaboration, 
Phys. Rev. D80, 054510 (2009), 0906.4720, 

[30] European Twisted Mass, P. Boucaud et al, (2008), 0803.0224, 

[31] European Twisted Mass, B. Blossier et al, JHEP 04, 020 (2008), 
0709.4574, 

[32] European Twisted Mass, C. Alexandrou et al., Phys. Rev. D78, 014509 
(2008), 0803.3190, 

[33] G. Martinelli, C. Pittori, C. T. Sachrajda, M. Testa, and A. Vladikas, 
Nucl. Phys. B445, 81 (1995), [hepdat/9411010| 

[34] P. Dimopoulos et al, PoS LAT2007, 241 (2007), 0710.0975, 

[35] P. Dimopoulos et al, Presentation at trento meeting on perspectives 
and challenges for full qcd lattoce qcd calculations, and in preparation, 
2008. 



[36] N. Mathur et al, Phys. Rev. D70, 074508 (2004), [Eep=pE70406 196 
[37] G. Meng and C. Liu, Phys. Rev. D78, 074506 (2008), 

10 



M. Luscher, Nucl. Phys. B364, 237 (1991), 

V. Bernard, M. Lage, U.-G. Meissner, and A. Rusetsky, JHEP 08, 024 
(2008), 0806.4495, 

SESAM, G. S. Bali, H. NefT, T. Duessel, T. Lipp ert, and K. Schilling, 
Phys. Rev. D71, 114513 (2005), [hepTat/0505012| 



P. Ball an d V. M. Braun, Phys. Rev. D54, 2182 (1996), [hep^ 
|ph/9602323[ 

A. P. Bakulev and S. V. Mikhailov, Eur. Phys. J. C19, 361 (2001), 



hep-ph/0006206, 



A. P. Bakulev and S. V. Mikhailov, Phys. Lett. B436, 351 (1998), 



hep-ph/9803298, 



K.-C. Yang, JHEP 10, 108 (2005), hep-ph/0509337 



J. A. Gracey, Phys. Lett. B488, 175 (2000), hep-ph/0007171 



K. G. Chetyrkin, J. H. Kuhn, and M. Steinhauser, Comput. Phys. 
Commun. 133, 43 (2000), [hep^ph/0004189| 

T. van Ritbergen, J. A. M. Vermaseren, and S. A. Larin, Phys. Lett. 
B400, 379 (1997), hep^ph/9701390 , 



M. Czakon, Nucl. Phys. B710, 485 (2005), [hep^piy0411261 



C. Hayne and N. Isgur, Phys. Rev. D25, 1944 (1982), 

C. Gattringer, L. Y. Glozman, C. B. Lang, D. Mohler, and S. Prelovsek, 
Phys. Rev. D78, 034501 (2008), 0802.2020, 

[51] G. Engel et al, (2009), 0910.2802, 



11 



